Improved Predictions of Reactor Antineutrino Spectra 



O 

(N 



X 
CD 

in: 



> 



X 



Th. A. Mueller/ D. Lhuillier/'B M. Fallot, ^ A. Letourneau/ S. Cornion,^ M. Fecliner,^ 
L. Giot,^ T. Lasserre,^ J. Martino,^ G. Mention,^ A. Porta, ^ and F. Yermia^ 

^Commissariat a I'Energie Atomique et aux Energies Alternatives, 
Centre de Saclay, IRFU/SPhN, 91191 Cif-sur-Yvette, France 
^ Laboratoire SUBATECH, Ecole des Mines de Nantes, Universite de Nantes, 
CNRS/IN2P3, 4 rue Alfred Kastler, U307 Nantes Cedex 3, France 
'"^Commissariat a I'Energie Atomique et aux Energies Alternatives, 
Centre de Saclay, IRFU/SPP, 91191 Cif-sur-Yvette, France 
(Dated: March 14, 2011) 

Precise predictions of the antineutrino spectra emitted by nuclear reactors is a key ingredient in 
measurements of reactor neutrino oscillations as well as of the recent applications to the surveillance 
of power plants in the context of non proliferation of nuclear weapons. We report new calculations 
including the latest information from nuclear databases and a detailed error budget. The first part of 
this work is the so-called ab initio approach where the total antineutrino spectrum is built from the 
sum of all /3-branches of all fission products predicted by an evolution code. Systematic effects and 
missing information in nuclear databases lead to final relative uncertainties in the 10 to 20% range. 
A prediction of the antineutrino spectrum associated with the fission of ^-^^U is given based on this ab 
initio method. For the dominant isotopes ^-^^U and ^"^^Pu, we developed a more accurate approach 
combining information from nuclear databases and reference electron spectra associated with the 
fission of ■^^^U, ^'^^Pu and ^*^Pu, measured at ILL in the 80's. We show how the anchor point of the 
measured total /3-spectra can be used to suppress the uncertainty in nuclear databases while taking 
advantage of all the information they contain. We provide new reference antineutrino spectra for 
235^^ 239p^ g^^^ 24ip^ isotopcs iu the 2-8 MeV range. While the shapes of the spectra and their 
uncertainties are comparable to that of the previous analysis of the ILL data, the normalization is 
shifted by about -1-3% on average. In the perspective of the re-analysis of past experiments and 
direct use of these results by upcoming oscillation experiments, we discuss the various sources of 
errors and their correlations as well as the corrections induced by off equilibrium effects. 



I. INTRODUCTION 

Nuclear power plants are the most intense man- 
controlled sources of neutrinos. With an average energy 
of about 200 MeV released per fission and 6 neutrinos 
produced along the /3-decay chain of the fission products, 
one expects some 2 x 10^° i^/s emitted in a Air solid angle 
from a 1 GW reactor (thermal power). Since unstable fis- 
sion products are neutron-rich nuclei all /3-decays are of 
/3~ type and the neutrino flux is actually pure electronic 
antineutrinos (i'e)- These unique features have been ex- 
ploited by several neutrino oscillations experiments 
Improvement in the accuracy of De spectra is motivated 
by next generation experiments aiming at unprece- 
dented sensitivity to the last unknown mixing angle ^13. 
The value of this parameter may determine the future 
trend of the neutrino physics, in particular for the search 
of CP violation in the lepton sector. Recent develop- 
ments of compact (1 m'^ target scale) S'e detectors as new 
safeguard tools for the monitoring of reactors would 
also benefit from an accurate description of De spectra. 

In a reactor core, only 1 neutron among the few gen- 
erated by the fission of a ^35^ nucleus should induce an- 
other fission, so that the core never reaches the over- 
critical regime. A fraction of the neutrons is actually 
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captured by the dominant ^ssy isotope leading to the 
production of new fissile isotopes: ^'^^Pu and to a lesser 
extent ■^'*^Pu. When operating, a core is thus burning 
^■^^U and accumulating ^^^^Pu. This is the so-called bur- 
nup process. In a pressurized water reactor, fission rates 
from both isotopes become comparable at the end of a 
cycle. The remaining fissions of ^*^Pu and fast neutron 
induced fissions of ^38^ share about 10% of the reactor 
power. As a result, the accurate prediction of the Ve spec- 
trum of a reactor requires following the time evolution of 
these four isotopes, as well as the knowledge of the asso- 
ciated /3-spectra of their neutron- rich fission products. 

This paper presents an improved treatment of the lat- 
ter piece of information, common to the prediction of 
the De spectrum of any moderated reactor. This work 
was triggered by the current effort of precision measure- 
ment of reactor neutrinos in the Double Chooz collab- 
oration 0]. Our approach combines the assets of the 
two main methods used so far. The first method is the 
so-called " ab initio approach" where the De spectrum as- 
sociated with one of the 4 fissioning isotopes is computed 
as the sum of the contributions from all fission products. 
This requires a huge amount of information on the thou- 
sands of /^-branches involved and the weighting factors 
of fission products, the fission yields. Section ITIl presents 
details on the ingredients of the ab initio approach while 
section HTll gives an update of ab initio calculations com- 
bining all data available today. The main systematic er- 
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rors are discussed and a prediction of ^ss^j gpgc^j-g, is 
given since this isotope is the only one with no integral 
beta spectrum measured yet. 

The second method relies on reference electron spec- 
tra j9l-[T2| measured at the high flux ILL reactor in Greno- 
ble (France) using a high resolution magnetic spectrom- 
eter [l^. It is presented in section ITVl where we explain 
how these electron spectra are converted into antineu- 
trino spectra with incomplete knowledge of the underly- 
ing physical distribution of /3-branches. We show how our 
" mixed-approach" can improve the control of systematic 
errors and lead to a significant correction of the refer- 
ence neutrino spectra used by all oscillation experiments 
so far. Finally we discuss in section IVl our results in the 
context of neutrino reactor experiments. 

II. INGREDIENTS OF REACTOR SPECTRA 

In the present work we describe the total j3 spectrum 
emitted by a reactor as the sum of the contributions from 
the four fissioning nuclei mentioned in section |T] 

Stot{E)^ akxSkiE) (1) 

where ak is the number of fissions of the fc"^ isotope at the 
considered time, Sk{E) is the corresponding /? spectrum 
normalized to one fission and E is the kinetic energy of 
emitted electrons. 

Most of the equations below can be found in textbooks 
but they are useful here to define our notations and dis- 
cuss the systematic errors in the following sections. In 
the ab initio approach, Sk{E) is broken up into the sum 
of contributions from all fission products. 

Sk{E) - J2 ^/f (^) (2) 

/p=i 

where Afp{t) is the activity of the fp*"^ fission product at 
time t and normalized to one fission of isotope " fc" . Then 
the spectrum Sfp{E) of each fission product is itself a 
sum of TVb /3-branches connecting the ground state (or 
in some cases an isomeric state) of the parent nucleus to 
different excited levels of the daughter nucleus 

Sfp{E) = ^ BR'j.p X Sfp{Zfp, Afp, E^fp, E) (3) 

6=1 

BR^fp and i^o/p ^^'^ branching ratio and the end- 
point energy of the 6*^ branch of the /p"^ fission product 
respectively. Zfp and Afp are the charge and atomic 
number of the parent nucleus. The sum of the branching 
ratios is normalized to the /3-decay partial width of the 
parent nucleus (1 if the parent is a pure /3~ emitter, < 1 
otherwise). 

Equations ([1]) to are valid for both electron and an- 
tineutrino spectra. The expression of the electron spec- 
trum of the 6*^ branch is given by the product of the 



following terms 

S% = X T{Zfp, Afp, E) X pE{E - E^fpf 

Norm. Fermi function Piiaso space 

X C'fpiE) X {l + S''fp{Zfp,Afp,E)) (4) 

Shape factor Correction 

To obtain the corresponding expression for the antineu- 
trino spectrum one can safely neglect the nucleus recoil, 
and replace in the above formula the electron energy E 
by the antineutrino energy 

E, = Elfp - E (5) 

By definition this one-to-one relation is valid only at the 
single /3-branch level. Thus this is a unique feature of 
the ab initio approach to predict electron and antineu- 
trino spectra with the same precision. The normalization 
factor K^jp of Eq.Q is calculated so that the integral 

/p^" Sfp{E)dE = 1. Hence the contribution to the in- 
tegral of Sfp{E) is driven by the branching ratio, as it 
should be. The next two terms come from the Fermi 
theory. The Fermi function J-{Zfp,Afp,E) corrects for 
the deceleration of the electron in the Coulomb field cre- 
ated by the Zfp x e positive charge of the parent nucleus. 
Therefore in the case of (3~ decay the Fermi function 
causes the electron spectrum to start at a non zero value 
at zero kinetic energy. This corresponds to a sharp step 
at the endpoint energy for the antineutrino spectrum, 
leading to discontinuities when summing several branches 
of different endpoints. 

The shape factor Cfp{E) brings extra energy depen- 
dence beyond the trivial phase space factor of the Fermi 
theory, due to the nuclear matrix element connecting the 
two nuclear levels of the /3-decay. Its complexity depends 
on the forbiddenncss of the transition, driven by the spin- 
parity of the connected levels. In the case of allowed 
transitions Cfp{E) is a constant and is absorbed in the 
normalization factor. 

For accurate predictions one must also take into ac- 
count corrections, represented by the (5^^ factor in Eq.®. 
This term is threefold 

S%iZfp,Afp,E) - Sqed{E) + Ac{Zfp,Afp) x E 

+ Awy-E (6) 

The 5qed term corrects for radiation of real and virtual 
photons by the charged fermion lines of the /3-decay ver- 
tex. Its expression has been calculated at order uqed by 
Sirlin et al. . The fact that only the charged fcrmions 
radiate photons implies that the Sqed formula differs for 
electron and antineutrino spectra, the electron spectrum 
deviating more from the shape predicted by lowest order 
calculation than that of the antineutrino. Strictly speak- 
ing, Eq.® now becomes Eq = Ee + E^ + E^ where Ej 
represents the energy of the radiated photon. Still the Ej 
spectrum goes like i/E^ and the dominant contribution 
comes from soft {E^ <^ Eq) radiated photons. Therefore 
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the total energy of the lepton pair remains very close to 
Eq . The physical constraint of conservation of the num- 
ber of particles is fulfilled by the equality 



SliE,) X (l + St^EDiE)) dE, 



' S)^{E,)y.[l + 6^QED{E))dE. (7) 



which we verified numerically. The Ac term is a Coulomb 
correction induced by the finite size of the decaying nu- 
cleus. It is related to the interference of < f? > and 
< a.r^ / > matrix elements, where R is the nuclear 
radius and a the spin operator. In the following we use 
the approximate expression derived by Vogel [l5| 
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and the Elton formula [l6[ as an estimate of nuclear radii. 

The nucleon itself also has a finite size and as a conse- 
quence the expression of its weak current deviates from 
that of a point-like particle. The complex effects of the 
nucleon internal structure can be absorbed in the defini- 
tion of factors in front of each term of the most general 
nucleon weak current allowed by the symmetries of the 
theory. The equivalent in the electromagnetic sector is 
the introduction of the Pauli (Fi) and Dirac (^2) form 
factors of the nucleon, with the F2 contribution propor- 
tional to the momentum transfer. The Aiy term con- 
tains the CVC partner contribution of F2 in the vector 
weak current and is called the weak magnetism correc- 
tion. Again we choose as a reference expression the one 
derived by Vogel [l5| 
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with A = gA/gv = —1.2695 the neutron disintegration 
constant and Mjq = 939 MeV the nucleon mass. One 
recognizes the (Kp — k„)/M7v term proportional to F2 
at small (^ -^at) square momentum transfer. Other 
expressions of Ac and A^/ can be found in the littera- 
ture [r^ . The associated uncertainty is large and ampli- 
fied by a sign compensation between the two terms. The 
net effect of these finite size corrections and its final error 
are discussed in sections UlI] and |lVl 



III. AB INITIO COMPUTATION OF BETA 
SPECTRA FROM FISSION FRAGMENTS 

A. Selection of the best data set 

In principle the ultimate prediction of the Sk{E) spec- 
tra comes from the knowledge of all quantities involved in 
equations to The first attempts at such ah initio 
approach were theoretical [l8l - l2l| . Most of these calcula- 
tions use rather crude models to describe the hundreds of 



involved nuclei but their goal is a correct description of 
total fission spectra Sk{E), not individual /3-transitions. 

Efforts have also been put recently in comparing micro- 
scopic models, mostly theoretical models based on QRPA 
and the nuclear Shell Model, in the framework of dou- 
ble beta decay studies [2^. Whereas these microscopic 
models are the ones susceptible to give the most reliable 
predictions, they are still difficult to apply to large sets 
of nuclei, especially heavy nuclei (such as the large mass 
region of the fission products) because of the large model 
spaces required. The estimation of the error associated 
to theoretical predictions remains a difficult task and in 
practice they are supplanted by measurements performed 
in the 80 's at the ILL High Flux Reactor in Grenoble [3- 
[T2j . Only the ^^^U spectrum remained calculated [20l.[2a| 
since no related data exist yet. Nevertheless a measure- 
ment in the fast neutron fiux of the FRMII reactor in 
Garching has lately been performed [23| and should be 
published soon. 

We describe here a complementary ab initio approach 
with the strategy of exploiting all data available in mod- 
ern nuclear databases while reducing the input of nuclear 
models. The total spectrum Sk{E) of each fissioning iso- 
tope is built up according to the equations of section HH 
retrieving the information on all /3-branches from the 
ENSDF nuclear database [Isj. The motivations for such 
an approach are that when all parameters of a /3-branch 
are known the neutrino branch is also known in a model- 
independent way and all errors on the input parameters 
can in principle be propagated. We have developed an in- 
terface with the ENSDF data library to read the relevant 
parameters of Eqs. ([3]) and ^ and their experimental 
error. The forbiddenness of a /^-transition is deduced 
from the spin and parity of the connected nuclear levels. 
In cases when this information is missing or uncomplete, 
the lowest possible forbiddenness is chosen by default. 
All transitions tagged as forbidden are then forced to 
be of unique type and the corresponding expressions of 
the shape factors C^p(i?) are polynomials in the electron 
and antineutrino momenta taken from [26l |. Using our 
homemade code BESTIOLE, all the above approxima- 
tions used for the calculation of each branch are tagged 
and various scenarios can be tested to estimate the error 
envelope of the final predicted spectrum. From Eqs.Q 
and ^ the electron and antineutrino spectra of each fis- 
sion product are computed and stored in a database. 

Then the total beta spectrum of one fissioning iso- 
tope is built as the sum of all fission fragment spectra 
weighted by their activity (Eql2]). These activities are de- 
termined using a simulation package called MCNP Utility 
for Reactor Evolution (MURE MURE is a precision 

code written in G-I--I- which automates the preparation 
and computation of successive MCNP (Monte-Carlo N- 
Particle transport code [l^) calculations either for preci- 
sion burn-up or thermal-hydraulics purpose. It is open- 
source, portable, and available at NEA [2^ and consti- 
tutes an efficient tool for non-proliferation and thermal 
power scenario studies (for more details see (30|). 
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FIG. 1. (Color online) Illustration of the detected antineu- 
trino spectrum in the case of ^^^U fissions (blue solid curve). 
Units are arbitrary and oscillation effects are suppressed. The 
detected rate rises from the threshold value at about 1.8 MeV, 
reaches a maximum around 4 MeV and vanishes after 8 MeV. 
This shape is the result of folding the emitted spectrum (black 
dashed-dotted curve), parameterization taken from section 
IV Al and beta-inverse cross section (red dashed curve). 



The detection process, common to many reactor an- 
tineutrino experiments, is the /3-inverse reaction on a free 
proton 



P 



(10) 



which sets an energy threshold for the antineutrino of 
1.804 MeV, the mass difference between the initial and fi- 
nal states (see figure [1]). Therefore the lowest energy part 
of the (3 spectra, below this threshold, is not addressed 
here. In this region, equilibrium is reached only after sev- 
eral months, requiring the control of significant transient 
effects when considering shorter irradiation times. Extra 
effects like the low energy /3-decays induced by neutron 



capture on ^■^'U and fission products [31[ would also have 



to be treated. On the high energy side, antineutrino rates 
above 8 MeV become negligible (< 0.5% of total detected 
rate) . This part of the spectrum is dominated by the very 
energetic (higli Q/s) transitions of rare exotic nuclei and 
cannot be accurately predicted. Thus the intermediate 
energy range resulting from the observation of the de- 
tected spectrum in figure [T] turns out to be favorable to 
the control of the systematic errors of the predictions of 
reactor antineutrino spectra. 

A powerful test of our calculations is the compari- 
son with the reference electron spectra from ILL [l0| - |12| . 
Such a consistency check gives valuable insight into the 
distribution of the numerous /3-branches, pointing to the 
main source of errors in the determination of the antineu- 
trino spectra. Considering all the data available in the 
ENSDF data library, the predicted /3-spectra associated 
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FIG. 2. Comparison of ^-^^U and '^'^^Pu reference electron 
spectra from 11] with our predictions based on the ab initio 
approach. The predictions have no free parameters and the 
rates are normalized to one fission. 



with the fission of ^ssy ^■^^Pu arc compared with the 
reference ILL data in figure [2l 

Although the spectrum falls quickly with energy, rea- 
sonable agreement is found on the shape and absolute 
normalization over a quite large energy range. Note that 
our prediction is parameter free. For finer analysis the 
residues of our predicted ^^Sy spectrum with respect to 
reference data are shown as the dashed-dotted line in fig- 
ured It reveals a ± 10% oscillation pattern of the calcu- 
lations around the data up to 7.5 MeV and a large overes- 
timation at higher energy. This overestimation points to 
the well known systematic effect of pandemonium [3^ . 
Indeed branching ratios and endpoints are usually de- 
termined by measuring the intensity and energy of 7- 
radiations emitted subsequently to the /? transition using 
high resolution but low efficiency Ge crystals. In the case 
of large Q/s a fraction of the beta branches connects the 
parent nucleus to very excited levels of the daughter nu- 
cleus. The strength of the associated low energy /3-rays is 
either spread over multiple weak 7-rays or concentrated 
in one high energy gamma ray. In both cases part or 
all the 7-cascade can be missed by the measurement ap- 
paratus. As a result low endpoint transitions are often 
missed and high endpoints are given too much weight in 
the global decay scheme of the parent nucleus. 

To correct for the pandemonium effect wc tried to 
gather /3-decay data using other experimental techniques 
than the /? — 7 coincidence. A valuable set of data comes 
from the measurement campaign undertaken by Teng- 
blad et al. [3^ in the late eighties at the on-line isotopes 
separators ISOLDE, at CERN, Geneva, and OSIRIS at 
the Neutron Research Laboratory, Studsvik. Some 111 
fission products, selected as the main contributors to the 
high energy part of reactor beta spectra (90% above 6 
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McV) were measured. Electron spectra were recorded 
independently from the emitted gamma rays. This pre- 
vented sensitivity to the pandemonium effect but at the 
same time part of the information on single /3-branches 
was lost. Among the 111 measured electron spectra, 44 
were found in perfect agreement with the spectra pre- 
dicted from the ENSDF database. The remaining 67 
were then replaced in our database. 
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FIG. 3. Residues of the -^^^U electron spectra computed as the 
difference of our ab initio calculations minus reference data 
from 11] divided by reference data. Dashed-dotted curve: 
ENSDF data only; dashed curve: some ENSDF data replaced 
by pandemonium corrected data; solid curve: unmeasured /3 
emitters are added on top of previous curve, using the gross- 
theory calculations of the JENDL nuclear database and few 
remaining exotic nuclei described by our model (see text). 

Another important source of data are the measure- 
ments based on Total Absorption Gamma Spectrometers 
(TAGS). The principle here is to implant the radioactive 
isotope on a foil surrounded by high efficiency gamma de- 
tectors able to collect the whole 7-cascade following the 
beta-decay. The distribution of total 7 energies gives ac- 
cess to the beta-strength of the studied isotope at the cost 
of a deconvolution analysis taking into account the full 
response of the apparatus. Eventually a complete beta- 
decay scheme can be determined providing the relevant 
beta-branch information for electron and antineutrino 



spectra. Thus the 29 nuclei of R. Greenwood et al. [34 1 
measured at the INEL facility, Idaho, were incorporated 
in our database. A. Algora and J. L. Tain studied care- 
fully both Tengblad et al. and Greenwood ct al.'s mea- 
surements [3^ . Both may be affected by several sources 
of systematic effects which are difficult to quantify. In 
particular both measurements of ^^Rb beta decay mean 
energy differ by more than 350 keV, while the ^^Rb decay 
scheme was used in Tengblad et al.'s analysis to quantify 
the 7-ray detector absolute efficiency. If the ^^Rb decay 
scheme is affected by the pandemonium effect, Tengblad 



et al.'s data sets may exhibit an overall systematic ef- 
fect. A new TAGS measurement of the decay properties 
of ^^Rb has recently been performed at the Jyvaskyla 
University facility |36j . and will help addressing the un- 
certainties associated to both sets of measurements. In 
cases when a fission product was present in both data 
sets (8 nuclei only), giving the priority to Greenwood et 
aZ.'s or Tengblad et aFs measurements changes the pre- 
dicted spectrum by 3% at most in the 4-5 MeV range, 
the effect drops at the 1% level or below elsewhere. In 
the following priority is arbitrarily given to Tengblad et 
aFs data. The dashed line in figure [3] shows the electron 
residues after merging our ENSDF based database with 
the above selected spectra supposedly corrected for the 
pandemonium effect. As expected the high energy part 
of our prediction has been significantly reduced leading 
to negative residues of increasing amplitude with energy. 
This indicates that a large part of the pandemonium ef- 
fect is probably corrected and that contributions from 
the missing unknown transitions of exotic nuclei grow 
rapidly with energy. 

To fill up this missing contribution we collected all 
available predictions of electron spectra from the JENDL 
nuclear database [13] . These predictions are based on the 
" Gross Theory of Beta Decay" fsH] and were included in 
the JENDL database to supplement ENSDF data show- 
ing incomplete level schemes or for nuclei for which no 
data were available. The estimated spectra were stored in 
the JENDL file so as to keep the consistency between the 
average decay energy value derived from the spectrum 
and that used for decay heat analysis [3§| . The calculated 
spectra were also compared with the directly measured 
spectra from the reference [s^ and revealed to be in very 
good agreement. The total contribution of the JENDL 
electron spectra, not already included in the ENSDF and 
pandemonium corrected data, was computed and con- 
verted to its associated total antineutrino spectrum fol- 
lowing the procedure described in |ll| . Finally the few 
remaining nuclei were described using a model based on 
fits of the distributions of the end-points and branching 
ratios in the ENSDF database, then extrapolated to the 
exotic nuclei. The result is the solid curve in figure [3] 
showing fiattened residues within a global envelope of ± 
10% over the whole energy range. This best agreement 
with the ILL reference data is actually valid for the ^'^^Pu 
isotope. Depending on the considered fissionning isotope, 
this new compilation of /3-decay data includes about 845 
nuclei and 10000 /3-branches, about 525 nuclei come from 
the ENSDF and pandemonium corrected data, 285 from 
the JENDL database and 35 from our model. It repre- 
sents the best data set for ab initio calculation. 



B. Results 

From this work we conclude that a compilation of all 
available data on the beta-decays of fission products can 
describe the antineutrino fission spectra at the 10% level. 
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illustrating the tremendous experimental work already 
achieved. Still, the relatively large energy range of de- 
tected antineutrino involves a sizeable contribution of 
unstable and poorly known nuclei in the total spectrum. 
Under these conditions, improving errors or even reach- 
ing the accuracy of the ILL reference spectra seems to 
require another fair amount of experimental effort. For 
applications like the determination of reactors decay heat 
calculations, a short list of " pandemonium candidates" to 
be remeasured with total absorption techniques has been 
identified Completion of a corrected beta-decay 

database is in progress (see for instance [4l[) with more 
and more refined analyses ji^]. Thanks to our database 
of fission product spectra, we have established a list of 
nuclei, contributing importantly to different energy bins 
of the antineutrino energy spectra from ^ssy g,^^ ^^gp^^ 
and that could be affected by the pandemonium effect. 
It will be the subject of a forthcoming publication. From 
this list we have selected a few nuclei which are amenable 
to experimental investigation using the TAGS technique, 
which can provide the beta intensity distribution in the 
full decay window eliminating the pandemonium effect. 
It appeared that some fission products being part of the 
measuremeiit priority list selected for reactor decay heat 
assessment [41| , also belong to the list of important con- 
tributors to the antineutrino emission in the energy win- 
dow of interest for neutrino oscillation studies. Recent 
and on-going experimental efforts carried out in the field 
of reactor physics, neutrino physics but also of interest 
for nuclear structure and astrophysics will certainly al- 
low to reduce the uncertainties associated to the reac- 
tor antineutrino spectra computed through the ab-initio 
method in the very next years [s^, IH, Hj] . We describe 
below the estimated error budget of our ab initio calcu- 
lations and give a prediction of electron and antineutrino 
spectra of ^^®U. In the perspective of neutrino oscillation 
analyses the ^35^ and ^^^Pu isotopes, which contribute 
to about 90% of a nuclear reactor spectrum, are predicted 
using a more accurate method presented in section IIVI 



C. Error Budget 

As mentioned earlier, the control of the parameters of 
all single /3-branches allows a full propagation of the er- 
rors quoted in the ENSDF database. All sources of error 
are treated as independent and the total error matrix of 
rates in energy bins is computed. In the simpler case 
of a spectrum at equilibrium, the activity of each fission 
product is approximated by the associated cumulative fis- 
sion yield indexed in the JEFFS. 1.1 database [i^. Then 
the uncertainty on branching ratios and fission yields can 
be propagated analytically while the uncertainty on end- 
points is propagated numerically (it turns out to have a 
negligible contribution). The dominant contribution of 
normalization errors induces large correlations between 
proximate bins as illustrated in table |T1 Note that these 
correlations are valid only for the specific part of the 
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TABLE L Correlation matrix in the range 2 to 3.5 MeV in 
250 keV bins obtained by propagating all sources of errors in 
ENSDF and JEFF databases. Branching-ratio errors cause a 
very high level of correlation then reduced by the end-point 
distribution and assumption of independent fission yields. 



measurement errors quoted in the nuclear databases. As 
summarized in table |TT] we know from the above section 
that systematic effects beyond these databases are dom- 
inant and will change these correlations in a non-trivial 
way as long as all /3-branches are not corrected. 



Kinetic 


Nuclear 


Forbid. 


Ac,w 


Missing 


E (MeV) 


databases 


treatment 


corrections 


info. 


2.00 


1.2 


0.2 


0.1 


10 


2.25 


1.3 


0.2 


0.2 


10 


2.50 


1.3 


0.1 


0.3 


10 


2.75 


1.3 


0.1 


0.3 


10 


3.00 


1.4 


0.4 


0.4 


10 


3.25 


1.6 


0.7 


0.5 


10 


3.50 


1.7 


0.1 


0.5 


10 


3.75 


1.9 


1.3 


0.6 


10 


4.00 


2.2 


1.6 


0.6 


10 


4.25 


2.5 


1.6 


0.7 


10 


4.50 


2.8 


1.4 


0.8 


10 


4.75 


3.2 


1.0 


0.8 


10 


5.00 


3.8 


0.5 


0.9 


10 


5.25 


4.4 


0.2 


0.9 


10 


5.50 


5.2 


0.2 


0.9 


15 


5.75 


6.1 


0.2 


0.9 


15 


6.00 


7.1 


0.2 


1.0 


15 


6.25 


8.0 


0.3 


1.0 


15 


6.50 


9.0 


0.4 


1.1 


15 


6.75 


10.1 


0.4 


1.1 


15 


7.00 


10.9 


0.5 


1.1 


20 


7.25 


11.0 


0.7 


1.1 


20 


7.50 


10.7 


0.8 


1.1 


> 20 


7.75 


11.1 


0.8 


1.2 


> 20 


8.00 


13.3 


1.2 


1.3 


> 20 



TABLE II. Sources of errors in the U electron spectrum as 
predicted by the ab initio approach. All errors are given in 
percent at la (68% CL). 



The second column of tabic |TI] lists the global effect 
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of the errors quoted in ENSDF at the 1 sigma level. It 
rises from 1 to 10 % in the 2-8 MeV range. Columns 
3 and 4 show the impact of the theoretical assumptions 
used to describe the shape of the /3-branches. Previous 
works always treated all branches as allowed. Comparing 
this hypothesis with our full treatment of forbiddenness 
shows changes of the final antineutrino spectrum below 
the 1% level (except for few bins around 4 MeV), vali- 
dating the allowed approximation. The error associated 
with the finite size corrections Ac and Aiy has been esti- 
mated by comparing the final spectra computed with no 
correction and those with the corrections from Vogel [lH| 
or Holstein p7j . The missing information on exotic nu- 
clei and the correction of the pandemonium effect unfor- 
tunately remain the dominant contribution in the final 
error of the ab initio approach. It is roughly estimated 
in the last column of the table based on the envelope of 
the various scenarios we tried and on the residues with 
respect to the reference ILL data. 




2 3 4 5 6 7 8 

Kinetic energy (IVIeV) 



FIG. 4. ab initio calculation of the electron (dashed his- 
togram) and antineutrino (solid histogram) spectra of ^'^^U 
for a 450 days irradiation time. For comparison the predic- 
tions of [l^] for an infinite irradiation time are plotted as 
crosses for the electron spectrum and stars for the antineu- 
trino spectrum. 



D. Predictions of U spectra 



Kinetic E 


Np 




t^ivie V ) 


(/fission/MeV) 


(/fission/MeV) 




12 h 


450 d 


12 h 


450 d 


z.UU 


1.13 


1.15 


1.43 


1.48 


Z.ZO 


9.81(-1) 


9.97(-l) 


1.26 


1.30 


z.oU 


8.47(-l) 


8.55(-l) 


1.12 


1.15 


Z. to 


7.24(-l) 


7.27(-l) 


9.80(-l) 


1.00 


O.UU 


6.11(-1) 


6.11(-1) 


8.70(-l) 


8.76(-l) 


o.zo 


5.07(-l) 


5.06(-l) 


7.57(-l) 


7.59(-l) 


O.OU 


4.16(-1) 


4.15(-1) 


6.40(-l) 


6.42(-l) 


O. ( 


3.37(-l) 


3.36(-l) 


5.39(-l) 


5.39(-l) 


4.UU 


2.68(-l) 


2.67(-l) 


4.50(-l) 


4.51(-1) 


4.Z0 


2.11(-1) 


2.10(-1) 


3.67(-l) 


3.67(-l) 


4.0U 


1.64(-1) 


1.63(-1) 


2.94(-l) 


2.93(-l) 


4. (0 


1.27(-1) 


1.27(-1) 


2.32(-l) 


2.32(-l) 


K nn 
O.UU 


9.72(-2) 


9.69(-2) 


1.83(-1) 


1.83(-1) 


o.zo 


7.37(-2) 


7.33(-2) 


1.43(-1) 


1.43(-1) 


cc m 


5.55(-2) 


5.52(-2) 


l.lO(-l) 


l.lO(-l) 


K ^K 
O.lO 


4.17(-2) 


4.14(-2) 


8.35(-2) 


8.35(-2) 


R nn 
D.UU 


3.12(-2) 


3.10(-2) 


6.21(-2) 


6.21(-2) 


D.ZO 


2.31(-2) 


2.30(-2) 


4.70(-2) 


4.70(-2) 


p. ccn 
D.OU 


1.68(-2) 


1.66(— 2) 


3.58(— 2) 


3.58(— 2) 


R 7K 


1.17(-2) 


1.16(-2) 


2.71(-2) 


2.71(-2) 


7.00 


7.92(-3) 


7.85(-2) 


1.95(-2) 


1.95(-2) 


7.25 


5.28(-3) 


5.23(-2) 


1.32(-2) 


1.33(-2) 


7.50 


3.48(-3) 


3.44(-2) 


8.65(-3) 


8.65(-3) 


7.75 


2.22(-3) 


2.19(-2) 


6.01(-3) 


6.01(-3) 


8.00 


1.40(-3) 


1.38(-2) 


3.84(-3) 


3.84(-3) 



TABLE III. ^-^^U electron and antineutrino spectra obtained 
by combining our best compilation of data sets with the activ- 
ity of all fission products as predicted by the MURE evolution 
code after a 12 h and 450 days irradiation time. Associated 
errors are those listed in table Ull 



applies to this prediction. Comparison with a previous 
estimate [l^l is illustrated in figure HI Despite some slight 
differences in shape, our work and previous predictions 
agree within ±10% across the full energy range. After 
multiplication with the /3- inverse cross section fEg fTO)) . 
the net effect on the integrated detected neutrino flux is 
a 9.8% increase for the ^ss^j contribution. 



As no experimental data on ^38^ g^j.g available at the 
present time we provide a prediction of its electron and 
antineutrino spectra using the above ab initio approach. 
This calculation has been performed using our best data 
set as defined in section llll Al The 238 u spectrum is given 
in table IIIII after two irradiation periods into a neutron 
flux: 12 h, similar to the irradiation time of the ILL data 
and 450 days as an approximation of a spectrum at equi- 
librium. All neutron captures effects are turned off for 
this prediction and the error budget described in table |TT] 



IV. IMPROVED CONVERSION OF REACTOR 
ELECTRON SPECTRA INTO ANTINEUTRINO 

In the previous section we showed that the ab initio 
approach had strong limitations due to unknown contri- 
bution from very unstable nuclei. Nevertheless we keep 
in mind that the beta-transitions described in nuclear 
databases represent about 90% of the total spectrum as 
measured at ILL. These physical distributions of end- 
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points and nuclear charges are precious information to 
control the conversion between electron and antineutrino 
spectra. We describe below how this can be combined 
to the very precise ILL electron data for an improved 
prediction of antineutrino spectra. 



A. Previous conversion procedure 

The measurements performed at ILL gave access only 
to the global electron spectrum of a fissile isotope, i.e. 
the sum of the contributions of all fission products. Thin 
target foils of fissile isotopes ^'^^U, ^'^^Pu and ^'^^Pu were 
exposed to the thermal neutron flux 80 centimeters away 
from the center of the compact fuel assembly. A tiny 
part of the emitted electrons could exit the reactor core 
through a straight vacuum pipe to be detected by the 
high resolution magnetic spectrometer BILL [l^. The 
electron rates were recorded by a point wise measurement 
of the spectrum in magnetic field steps of 50 keV, provid- 
ing an excellent determination of the shape of the electron 
spectrum with sub-percent systematic error. The pub- 
lished data are smoothed over 250 keV. Except for the 
highest energy bins with poor statistics, the dominant 
error was the absolute normalization, quoted around 3% 
(90% CL) with weak energy dependence. Note also that 
the ILL spectra are taken after typically 1 day of irradia- 
tion, meaning that the longest lived beta emitters (lowest 
energy beta rays) haven't reached equilibrium yet. These 
aspects are discussed in detail in section IVB] 
The neutrino spectra, not directly detected, were de- 
duced from those of the electron via a conversion proce- 
dure which induced some extra systematic effects. In 0- 
[T^ the authors considered 30 virtual beta branches. The 
procedure consisted in dividing the electron spectrum 
into 30 slices. Starting with the highest energy slice, the 
few data points in this slice were used to fit the endpoint 
and branching ratio of the first virtual branch. The full 
contribution of this virtual branch (from endpoint down 
to zero energy) was then subtracted from the experimen- 
tal spectrum and the procedure repeated for the next, 
lower energy, slice. Then the antineutrino spectrum was 
simply the sum of all fitted virtual branches, converted to 
antineutrino branches by replacing by E^, ~ Eq — E^. 
and applying the correct radiative corrections. This pro- 
cedure was repeated several times, describing the spec- 
trum with somewhat difl[erent sets of end point energies. 
Possible steps in the shape induced by the relatively small 
number of virtual branches were smoothed out by taking 
the average of all spectra and merging the 50 keV bins 
into the 250 keV presented in the publications. The the- 
oretical expression of a virtual branch was the same as 
Eq. (|4]) , except for the Ac and Aw corrections which were 
treated at the very end as an effective linear correction 
to the final antineutrino spectra 

AN^^{E,) ~ 0.65 {E, - 4.00) % (11) 



with El, in MeV. The final error of the conversion pro- 
cedure was estimated to be 3-4% (90% CL), to be added 
in quadrature with the electron calibration error. 

In a recent paper [46l |. P. Vogel pointed out the main 
limitations of this conversion procedure. Despite the dis- 
continuity of the Fermi function at the endpoint energy of 
an antineutrino branch, the " true" antineutrino spectrum 
from fission fragment appears continuous because thou- 
sands of branches contributes with a quasi-continuous 
endpoint distribution. When describing the spectrum by 
only 30 virtual branches a spurious oscillation with re- 
spect to the true spectrum appears around each virtual 
endpoint energy. Therefore smoothing out these oscilla- 
tions requires sufficiently narrow slices of electron data 
and antineutrino energy bins several times larger than 
the slice width. All these criteria couldn't be fulfilled for 
the electron data taken at ILL and the estimation of the 
remaining effects is part of the quoted error bar. The 
other criterion highlighted by P. Vogel was the knowl- 
edge of the average nuclear charge < Z > of the virtual 
branches as a function of their endpoint energy. This in- 
formation turns out to be of crucial importance for the 
shape of the high energy part of the antineutrino spec- 
trum. 



B. Improvements of the conversion procedure 

Our new conversion method allows us to address the 
sources of errors in a complementary way. It consists 
in starting with our ab initio prediction of section IIIII 
and restrict the use of effective branches to fit only the 
missing few percent contribution of the difference with 
the reference ILL electron data. This way we keep the 
distributions of beta branches very close to the physi- 
cal one and we can apply Ac and A]y corrections at the 
branch level. The reference ILL electron data are still fit- 
ted but the contribution of unphysical virtual branches 
is reduced by an order of magnitude. We use all avail- 
able data in ENSDF plus the above-mentioned 67 nuclei 
from the "pandemonium corrected" measurements. In 
the case of the data from reference [s^l only the total /3 
spectra of each nuclide are available, not the complete 
decay scheme as would be best. Hence, to be converted 
to a neutrino spectrum, each /3 spectrum measured must 
be fitted by a set of branches. These branches differ from 
the virtual branches used to fit the ILL data by the fact 
that in this work the nuclear charge is perfectly known. 
Moreover the validation of the spectrum shape at the 
level of one nucleus provides a more refined description 
of the /3-decay scheme equivalent to a small slice size in 
the total spectrum, pinning down the sources of errors 
discussed in f46j . 

On top of the contributions of all ENSDF and pande- 
monium corrected /3-branches, the missing contribution 
to match the ILL electron spectrum is fitted using a set 
of 5 effective /3-branches with a nuclear charge of Z=46 
(chosen as the average of the distribution of fission prod- 
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FIG. 5. (Color online) The blue hatched area shows the con- 
tribution of our ab inito prediction (ENSDF + pandemonium 
corrected nuclei) relative to the ILL reference data. The miss- 
ing contribution coming from unknown nuclei and remaining 
systematic effects of nuclear databases (red hatched area) is 
fitted using a set of 5 effective /3-branches. 



ucts), and assuming that transitions are allowed tran- 
sitions. The normalization and the end-point are two 
free parameters for each branch. An example of the " ab 
initio'^ and "fitted" contributions of ^ssy electron spec- 
trum is illustrated in figure [5] as stacked histograms. The 
residues of the fitted missing contribution are shown in 
figure [6ja). They are small, typically at the level of the 
statistical error of the ILL data, except in the 4.5-6.0 
MeV range where one can see an oscillation pattern with 
an amplitude reaching three times the error of the ILL 
data at maximum. This may point to a systematic ef- 
fect due to a failure of the fit model. We checked that 
using more effective branches is not efficient because of 
the limited number of experimental points available. The 
impact of these non statistical residues in the final error 
is discussed latter. 



C. New reference antineutrino spectra 

Converting all branches from the nuclear databases 
plus the 5 fitted ones into antineutrino branches (as de- 
scribed in section |ll| , we obtain the predicted antineu- 
trino spectrum. The residues with respect to the pre- 
diction of Schrcckenbach et al. arc shown in figure [ni^b). 
It exhibits a good agreement in shape but a mean nor- 
malization shift of about -1-3%. This shift of the emit- 
ted antineutrino flux is modulated at higher energy by 
oscillations which look like images of the oscillations in 
the electron residues. When folded with the /3-inverse 
cross section the predicted increase of detected antineu- 
trino rate from ^35^ is about 2.5%. Note that this posi- 
tive shift is a mean value computed between the energy 
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FIG. 6. Relative difference to ILL spectra for electron (a) and 
antineutrino (b) spectra for ^^^U. While electron residues are 
zeroed at the ±1% level by the fitting procedure, antineu- 
trino residues exhibit a mean normalization shift of about 3% 
(dashed line in(b)). 



threshold (1.8 MeV) imposed by the detection process 
and infinity. The physical constraint of one emitted elec- 
tron for one emitted antineutrino must still be fulfilled. 
We did check that the total integral of our electron spec- 
trum fitted on the ILL reference and the integral of the 
associated converted neutrino spectrum were identical at 
the 10"'' level. Since all our individual (3 and antineu- 
trino branches are normalized to the same integral with 
much higher accuracy this result gives an estimate of the 
numerical precision in the sum of thousands of branches. 

To test the validity of our procedure we applied it to ef- 
fective calculated electron and antineutrino spectra. This 
method was inspired from the work of P. Vogel [1^ . We 
generated electron and antineutrino spectra as the sum 
of the spectra of all fission products indexed in ENSDF, 
weighted by the activity predicted by the MURE code 
after 12 hours of irradiation. We know from section Hill 
that these spectra are close to the ones measured at ILL 
and in the context of this test we call them "true" spec- 
tra in the sense that they are unambiguously connected 
to each other by the conversion of each single branch of 
the sum. Then we followed the exact same procedure as 
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FIG. 7. (Color online) Independent cross-check of our results 
based on known reference spectra (pure ENSDF database). 
Dashed red line: electron residues after fitting with 30 vir- 
tual branches. Dotted blue line: relative difference between 
the reference antineutrino spectrum and the one converted 
according to the ILL procedure, in 50 keV bins. Smoothing 
out the residual oscillations in 250 keV bins (solid black line) 
exhibits the same 3% normalization shift than in figure [6jb). 



the one described in |10l - ll2l | to convert our "true" elec- 
tron spectrum into an antineutrino spectrum. This in- 
cludes using 30 virtual branches and the same effective 
Z distribution and the same effective Ac + Ay/ correc- 
tion (Eg.pTjl). The spectrum converted in this fashion 
is finally compared to our " true" antineutrino spectrum. 
Figure [7] shows that despite a very good fit quality of the 
electron spectrum (all electron residues are within a few 
10"'^ from 1 to 8 MeV) the converted antineutrino spec- 
trum exhibits residues with oscillations of few percent 
amplitude around the endpoint of each fitted branch. As 
expected, rebinning smoothes out these oscillations (solid 
curve) but a residual ~ -1-3% offset is clearly visible across 
the whole energy range. This curve can be directly com- 
pared to the result of our conversion of the ILL data, 
in figure [6l[b) . Very good agreement is found, validating 
the above predicted deviation from the ILL antineutrino 
spectra. 

Switching on and off the various ingredients of the ILL 
conversion reveals a twofold origin of the normalization 
shift. At low energy, the deviation is mainly due to the 
treatment of the Ac and Ay/ correction terms. In Eq.([6|) 
these two terms appear multiplied by the energy hence 
one would expect their net contribution to grow with 
energy. Nevertheless Ac itself has some hidden energy 
dependence via the Z distribution of all branches (EqE]) 
while the estimate of Ay/ is a constant (Eql9]). To illus- 
trate the size of the total correction Ac + Ay/ we show 
the ratio of corrected over non-corrected spectra in figure 
[S] in the case of total spectra build up from all ENSDF 
branches. There is a linear trend at high energy but di- 
rect comparison with the effective linear correction used 
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FIG. 8. Ratio of total spectra computed from all ENSDF 
branches with and without the Ac and Aw correction terms. 
The dashed (solid) curve shows the ratio of electron (neutrino) 
spectra. 



in previous analysis Eq. (fTT|) shouldn't be made at this 
stage. In fact in the conversion procedure the electron 
spectrum is fitted on the electron ILL data. By defini- 
tion the fitting procedure optimizes the parameters of the 
few virtual electron branches used in this work so that 
the total electron spectrum always matches the ILL data, 
whatever the correction terms included in the theoretical 
expression of the branches. Therefore in the conversion 
process only the neutrino spectrum is sensitive to the Ac 
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FIG. 9. The solid curve is the same solid curve as in figure [T] 
The dashed curve illustrates the deviation induced by the 
implementation of Ac and Aw corrections at the level of each 
virtual branch. It explains most of the -1-3% shift observed in 
figures [Bl^b) and [7] below 4 MeV. 
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and Ayy terms. The final effect is not intuitive but ean 
be computed numerically as shown in figure [9l When 
applied at the branch level the Ac and Aw corrections 
deviate from the effective formula of Eq. (fTTj) and explain 
the predicted 3% shift in the low energy region. We stress 
that no improvement was brought to the theoretical ex- 
pression of the corrections (Eq. [5]and[S]), only the way of 
implementing these corrections has been revisited avoid- 
ing the extra approximation of the effective correction 
of Eq. pTj) . While justified at high energy, we show that 
this approximation was giving too much amplitude to the 
correction at low energy. The effect is large enough to 
explain the observed 3% shift at low energy. 
At higher energy the dominant source of the shift comes 
from the parameterization of the charge distribution as- 
sociated with the virtual /3-branches. Figure [TU] illus- 
trates how large errors can be induced by the too rough 
approximation of a constant nuclear charge. In the ILL 
data analysis the mean charge of each virtual branch was 
taken from a polynomial fit Z(Eq) of the tabulated nu- 
clear data (Eq (3) of [l^). This greatly reduced the 
bias but still this approach doesn't take into account the 
very large dispersion of nuclear charges around this mean. 
Even a new Z{Eq) function, fitted on the same data used 
to build our "true" spectrum generates a deviation of a 
few percent at high energy. 

Further cross-checks of our results are based on min- 
imizing the electron residues in different independent 
ways. First we checked that the new conversion proce- 
dure is not sensitive to the chosen starting point for the ah 
initio calculations. All the results presented above were 
obtained by adding 5 virtual branches to the spectrum 
built up from nuclear databases and using independent 
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FIG. 10. Deviation from the true antineutrino spectrum in- 
duced by various Z(E) functions used in the formula of the 
virtual branches. For this test Ac and Aw are turned off in 
the "true" branches as well as in the virtual branches. This 
effect takes over the +3% deviation observed in figures [Bj^b) 
and □ above 4 MeV. 
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FIG. 11. Stability check of antineutrino residues when using 
different estimates of fission yields for the initial ab initio cal- 
culation (see eq.[2]): independent yields calculated via MURE 
after 12h of irradiation (solid line), after 36h (dotted line) and 
cumulative yields (dashed line). The difference between neu- 
trino residues is less than 1% below 4.5 MeV and negligible 
for higher energies 



yields calculated after 12 hours irradiation time to match 
as closely as possible the experimental conditions at ILL. 
Figure [TT] shows the variation of the neutrino residues 
when using independent yields at 36 hours instead of 12 
hours or even cumulative fission yields, corresponding to 
the equilibrium regime reached after infinite irradiation 
time. One can see that the variations induced in the an- 
tineutrino spectra are negligible (< 1%). This can be 
understood by the fact that although the different sets of 
fission yields do change the shape of the ah initio spec- 
trum by a few percent, this change is absorbed by the 
virtual branches fitting the missing contribution with re- 
spect to the ILL electron spectrum. The underlying dis- 
tributions of nuclear charges and end-points remain very 
similar leading to the same final residues. This illustrates 
how our mixed approach gets rid of the dominant errors 
of the ah initio approach. 

In order to avoid the use of virtual branches we also 
tried minimizing the electron residues by tweaking the in- 
put parameters of the ah initio calculation, namely, the 
distributions of branching ratios or end-points or both at 
the same time. This technique makes sense only if the 
tweaking docs not disturb too much the physical distri- 
butions. Therefore in the minimization procedure we im- 
plemented limitations of typically 15% for the variation 
range of physical parameters. Electron residues of similar 
quality could be obtained in this way at low energy but 
deviations of several percents couldn't be avoided above 5 
MeV. As already observed in figure [HI any large residues 
pattern in the electron fit shows up in the antineutrino 
residues, slightly shifted in kinetic energy and amplified 
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TABLE IV. Results of the new conversion procedure on the 
^^^U antineutrino spectrum. The electron residues between 
our prediction and ILL data are given in percent as an indica- 
tion of the quality of the fitting procedure. The antineutrino 
spectrum is normalized per fission and correspond to the spec- 
trum for a 12h irradiation time. All errors are given in percent 
at la (68% CL). 



FT 

1 ivie V 1 


p res 


239p^ 

(/fission) 


A AT 
[/o) 


p res 
[/o) 


AT 

i^l fission) 


A AT 


o nn 
z.UU 


n nQ 
-U.Uo 


1. 


13 


Q 
Z.O 


n c\A 
-U.U4 


1 07 
i.Z / 


O 

z.z 


o oc; 

Z.ZO 


n 1 o 
-U.iz 


9.19( 


-1) 


O Q 
Z.O 


n 1 n 
U.iU 


1 n7 
i.U / 


o o 
Z.Z 


o p;n 
z.OU 


C\ ^ A 
U. 14 


7.28( 


-1) 


A 
Z.4 


nil 


y.uo(— i ) 


o 
Z.Z 


O IK 

z. /O 


n QQ 
-U.OO 


6.13( 


-1) 


A 
Z.4 


n nn 
U.UU 


7 fiQ / 1 ^ 


o 
Z.Z 


Q nn 
o.UU 


n Qi 
U.oi 


5.04( 


-1) 


O A 

Z.4 


n OA 
U.Z4 


Qn/ 1 ^ 

o.oy[— i ) 


o o 
Z.Z 


Q op; 
o.zO 


n np; 
U.UO 


4.10( 


-1) 


A 
Z.4 


n nQ 

u.yo 


p; QI / 1 ^ 


o 
Z.Z 


Q p;n 
O.OU 


n c\A 
U.U4 


3.21( 


-1) 


a 

Z.O 


n fiK 
U.oO 


/I QQ / 1 ^ 

4.ooi^ — i ) 


/I 
z.4 


6.(0 


i.4y 


2.54( 


-1) 


a 
Z.O 


n AG 

u.4y 


Q p;i / 1 ^ 


/I 
z.4 


A nn 
4.UU 


n QT 
-U.o / 


2.00( 


-1) 


O T 

Z. / 


n no 

-U.uy 


o Qo/ 1 ^ 
Z.oZl^ — i ) 


o p; 
Z.O 


A op; 
4.zO 


n fiQ 
-U.UO 


1.51( 


-1) 


n 
z.y 


n no 
U.UZ 


o 1 Q / 1 ^ 


7 
Z. / 


A p;n 

4.0U 


4.4y 


1.10( 


-1) 


Q n 
o.u 


1 Ofi 
-l.ZO 


1.001^ — i ) 


Q 
Z.O 


A IK 
4. /O 


n 01 
U.zi 


7.97( 


-2) 


Q n 
o.u 


n Qn 
-U.oU 


1 00/ 1 ^ 
i.zzy — i ) 


Q 
Z.O 


p; nn 
o.UU 


Al 
-Z.4 / 


6.15( 


-2) 


Q Q 

o.o 


n /IS 

-U.4o 


Q p;q/ 0"^ 
y.oy — z ) 


Q 1 
0. 1 


p; op; 
O.zO 


O /I Q 
-z.4o 


4.68( 


-2) 


Q Q 
O.O 


n no 
-u.yz 


7 Q(^/ o^ 
/ .oO( — Z ) 


Q 1 
O.l 


p; p;n 
O.OU 


K A^ 


3.50( 


-2) 


A A 
4.4 


n nQ 
-U.yo 


K K0( 0^ 
O.OZl^ — Z ) 


A Q 
4.0 


p; IK 
0. (0 


n Qo 
-U.oz 


2.55( 


-2) 


A fi 
4.0 


n nv 
-U.U / 


/I ni / o^ 


A K 
4.0 


(\ nn 
D.UU 


O 0^^ 
Z.ZO 


1.82( 


-2) 


A n 
4.y 


1 «G 

i.oy 


o QI / o^ 

Z.oil^ — Z J 


A 7 
4. / 


f\ OK 
D.ZO 


1 p;q 
i.oo 


1 QO ^ 

i.oZi 


-2) 


p; n 
o.u 


n 77 

U. / / 


on/1/ 0^ 

Z.U4l^ — Z J 


A 7 
4. / 


C\ KC\ 
D.OU 


n IP-. 
-U. /O 


9.82( 


-3) 


p; o 
O.Z 


n 1 n 
U.IU 


1 p;n/ o^ 


A n 
4.y 


IK 
D. / 


K Al 
0.4 / 


7.32( 


-3) 


p; o 
O.Z 


n n/1 
-u.y4 


1 n7/ o^ 
i.U / 1^ — Z J 


A O 

4.y 


7.00 


-5.01 


5.13( 


-3) 


7.1 


-0.32 


7.20(-3) 


5.3 


7.25 


-1.73 


3.15( 


-3) 


9.2 


-1.23 


4.47(-3) 


5.3 


7.50 


-8.94 


1.83( 


-3) 


11.1 


-0.96 


2.54(-3) 


5.7 


7.75 


-32.75 


1.03( 


-3) 


15.7 


-1.07 


1.65(-3) 


5.7 


8.00 


-55.56 


4.91( 


-4) 


20.6 


-1.55 


9.63(-4) 


7.0 



TABLE V. Results of the new conversion procedure on Pu 
(1.5 days irradiation time) and ^*^Pu (1.8 days irradiation 
time) antineutrino spectra - see comments for ^^^U. The con- 
version error and the error due to the Ac,w terms are the same 
as those given in table HVl The statistical and normalization 
errors for both plutonium isotopes can be found in 12]. 



by a factor of about 3. This behavior is used to estimate 
our error budget, summarized in table IIVI Even in the 
case of zero residues, we count the statistical error of the 
reference /3 spectrum as the minimum error of the con- 
verted antineutrino spectrum (column 4). Based on our 
numerous tests of fitting methods, the amplified envelope 
of non statistical electron residues is added in quadrature 
as an extra conversion error (column 5). 

In column 6 the error due to the Ac,w terms is 
computed by propagating a 100% relative uncertainty 
through the conversion procedure. Finally the normal- 
ization error of the ILL reference data is directly propa- 
gated as a normalization error of the converted antineu- 
trino spectrum (column 7). The total error is taken as 
the quadratic sum of all previous sources of errors. In the 
perspective of neutrino oscillation analyses it is manda- 
tory to consider the correlations between energy bins. 



The statistical and conversion errors are driven by ran- 
dom processes. Therefore they do not induce any bin 
to bin correlation. The normalization error of the ILL 
data should be treated as fully correlated over the whole 
energy range. Regarding the Ac^w terms we observe 
that they are propagated as a linear correction to the 
converted antineutrino spectrum above > 4 MeV. The 
uncertainty on the slope coefficient fully correlates all 
high energy bins. Below 4 MeV the precise determina- 
tion of the correlations would require dedicated numerical 
studies but in that energy domain the size of the Ac^w 
corrections and their error are small and have negligible 
impact in the error budget. 

We applied our same conversion procedure to pluto- 
nium isotopes. Results are given in table |Vl The con- 
clusions given for the 235]j antineutrino spectrum remain 
valid for these isotopes. The main net effect is a mean 
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« +3% normalization shift with respect to previous ref- 
erence spectra. The equivalent increase in the detected 
neutrino spectrum is 3.1% for ^'^^Pu and 3.7% for ^^^Pu. 

V. CONTEXT OF REACTOR NEUTRINO 
EXPERIMENTS 

A. A useful phenomenological parameterization 

A phenomenological parameterization of our fission an- 
tineutrino spectra could be useful for sensitivity studies 
requiring different binning or energy domains than those 
proposed in tables [nil lY] and El Therefore, as in [13 
we provide a parameterization of the spectrum of a given 
isotope using the exponential of a polynomial 

6 

SkME.) = exp (^apfe^r^). (12) 

p=i 

with the coefficients apk determined by a fit to the data 
using the MIGRAD algorithm of the TMinuit ROOT 
class fZsj . To this aim we minimize the x^-function 

hi 

6 

with D,=Y1 (^pkiEl'^f - In^r' (13) 
p=i 

where eI"^ and sl!^ = Skisi'^) are the values of the 
antineutrino energy and the corresponding antincutrino 
spectrum, respectively, provided in tables IIIIl HVl and El 
Since we are fitting the logarithm of the flux the co- 
variance matrix Vij contains the relative errors of the 

Sj, . For the diagonal element Vu we take the total error 
quoted in the above mentioned tables. Because the error 
of the 238^ spectrum has been estimated from the en- 
velope of all systematic effects of the nuclear databases, 
we assume no bin to bin correlation. For our new con- 
verted spectra, the fully correlated errors on the absolute 
calibration of ILL /? spectra and on the Ac.w correction 
terms contribute to the off-diagonal elements of the co- 
variance matrix as 

T / ^cal cal I ^corr ^corr ^ -L A / 1 /I ^ 

Vij = cr,; aj aj , i j. (14) 

The plots (a) to (d) of figure [T2l shows the resulting spec- 
tra for the 6 parameter fits in comparison to the data and 
the corresponding x^-values per degree of freedom. The 
best goodness-of-fit is obtained with polynomial of order 
five. The plots (e) to (h) of figure [12] show the residues of 
the fit in units of Ui , where the error is obtained from the 
covariance matrix V by Oi = S^. WVu- Note that in case 
of correlations between the these residuals do not add 
up to the total x^- All antineutrino spectra are very well 
described by the chosen phenomenological parameteri- 
zation of Ea. (fT^ . The best fit coefficients apk and their 



correlation matrix are given in table IVTl We can see quite 
large anticorrelations among consecutives fit parameters 
which could be induced by the choice of the exponential 
of a polynomial for the fit function. Therefore one should 
be aware of possible bias in the propagation of correla- 
tions when using these fits whereas it is a practical way 
to compute nominal spectra. 
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TABLE VI. Coefficients apk of the polynomial of order 5 for 
antineutrino flux from elements k = ^^Sy, ^ssy^ 239p^ ^^^^ 
2^^Pu. In the column Sapk the la error on apk are given. 
Furthermore the correlation matrix of the errors is shown. 
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FIG. 12. (Color on line) Illustration of the fit of the antineutrino spectra predicted in this work. The red curves in the left 
panels correspond to a 6 parameters fit (polynomial of order 5). Also shown the data with their la error bars and the per 
degree of freedom. In the right panels we show the residual of the fits. 
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B. Off-equilibrium corrections 

The ILL spectra were acquired after a quite short ir- 
radiation time in a quasi pure thermal neutron flux, be- 
tween 12 hours and 1.8 days depending on the measured 
isotopes. For neutrino reactor experiments the irradia- 
tion time scale would rather be a reactor cycle duration, 
typically 1 year. Among the fission products, about 10% 
of them have a /3-decay life-time long enough to keep ac- 
cumulating after several days some of them presenting 
sufficient large capture cross sections to possibly affect 
the final inventory. Moreover, in a standard PWR, the 
neutron energy spectrum exhibits more important ep- 
ithermal and fast neutron energy components than in 
the ILL measurements. These higher energy components 
of the neutron flux add small contributions to the fis- 
sions of 235u (|g^g ^f^i\ f-Qj. i^jj^Q other fissioning isotopes) 
leading to different distributions of the fission products. 
In this section we study the effect of these phenomena 
on the reference neutrino spectra and compute the as- 
sociated corrections. Since these corrections are relative 
deviations between spectra at different irradiation times, 
we assume they are pretty insensitive to the sources of 
error of our ab initio calculations discussed in section Hill 

Therefore the study was done with the MURE simu- 
lation of a PWR assembly of N4 type exhibiting a mod- 
eration ratio equal to the one of a PWR core in order 
to represent the full reactor neutronic conditions. The 
infinite multiplication coefficient of the simulation has 
been successfully compared with similar simulations per- 
formed with the deterministic code DRAGON for french 
PWRs [i^ . This simulation represents thus very well the 
real physical conditions of a reactor core. Condition of 
a constant power is assumed, renormalizing the neutron 
flux at each time step in order to compensate for the 
fuel burnup. We adapted the code in order to compute 
and store the amount of all /3~ emitters produced over 
time jsoj . In our simulation, the fission yields from the 
JEFFS. 1.1 nuclear data library [i^ were used. The yields 
coming from the 25 meV, 400 keV and 14 MeV libraries 
were weighted by the fission rates in each neutron energy 
region. This simulation was compared with independent 
calculations using the FISPACT code [5l| based on the 
EAF nuclear data library which just evolves the isotopic 
concentrations over time. A constant mean neutron flux 
of 3.10"'^"'neutron/cm^/s was used for the FISPACT cal- 
culations. 

The departures from our reference spectra are dis- 
played as a function of time in table IVIII for some rel- 
evant low energy bins. As expected, the accumulation of 
long-lived nuclei shows up as positive deviations which 
amplitude decreases with the neutrino energy and be- 
comes negligible above 3.5 MeV. At the threshold of the 
beta-inverse process it takes about 100 days of irradia- 
tions for the antineutrino spectrum to be stable at the 
1% level. Noting that the irradiation time for the refer- 
ence spectrum of ^asy jg 12 h instead of 36 h for ^'^^Pu 
and ^^^Pu, the corrections are similar for all isotopes. 



235 ^ 





2.0 MeV 


2.5 MeV 


3.0 MeV 


3.5 MeV 


4.0 MeV 


36 h 


3.1 


2.2 


0.8 


0.6 


0.1 


100 d 


4.5 


3.2 


1.1 


0.7 


0.1 


1E7 s 


4.6 


3.3 


1.1 


0.7 


0.1 


300 d 


5.3 


4.0 


1.3 


0.7 


0.1 


450 d 


5.7 


4.4 


1.5 


0.7 


0.1 






2.0 MeV 


2.5 MeV 


3.0 MeV 


3.5 MeV 


4.0 MeV 


100 d 


1.2 


0.7 


0.2 


< 0.1 


< 0.1 


1E7 s 


1.3 


0.7 


0.2 


< 0.1 


< 0.1 


300 d 


1.8 


1.4 


0.4 


< 0.1 


< 0.1 


450 d 


2.1 


1.7 


0.5 


< 0.1 


< 0.1 


241p^ 




2.0 MeV 


2.5 MeV 


3.0 MeV 


3.5 MeV 


4.0 MeV 


100 d 


1.0 


0.5 


0.2 


< 0.1 


< 0.1 


1E7 s 


1.0 


0.6 


0.3 


< 0.1 


< 0.1 


300 d 


1.6 


1.1 


0.4 


< 0.1 


< 0.1 


450 d 


1.9 


1.5 


0.5 


< 0.1 


< 0.1 



TABLE VII. Relative off-equilibrium correction (in %) to 
be applied to the reference antineutrino spectra listed in ta- 
bles |IV| and |Vl for several energy bins and several irradiation 
times significantly longer than the reference times (12h U for 
and 36h for Pu). Effect of neutron captures on fission prod- 
ucts are included and computed using the simulation of a 
PWR fuel assembly with the MURE code. 



We checked with our evolution codes that the effects of 
neutron capture on the fission products as well as the 
contribution of the neutron spectrum above the thermal 
energy domain have small impact on the off-equilibrium 
corrections. Other tests have been performed with the 
FISPACT code, showing that these results may depend 
on the system (neutron, flux and energy spectrum, geom- 
etry) used in the calculation. The error envelop covering 
our different results is of 30% on the total off-equilibrium 
corrections. Therefore the results quoted in table IVIII 
should be taken as typical corrections at a N4 reactor. 
For applications with signicantly different neutron flux 
or fuel geometry, dedicated simulations should be car- 
ried out for an accurate correction of the lowest energy 
bins of the antineutrino spectrum. 

Off-equilibrium effects have independently been evalu- 
ated for the analysis of the Chooz experiment [13 , which 
measured the neutrino spectrum of the two N4 reactors 
of the Chooz site. In this reference, the departure from 
the antineutrino ILL spectra were computed using the 
cumulative yields of some known long-lived fission frag- 
ments. The results are shown as markers in figure [13] 
to be compared with the histograms of our calculations. 
The overall agreement is good, even when evolving the 
spectrum back to irradiation time as short as 10^ s, where 
the corrections become quite large and have steep varia- 
tions in time. 
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atic errors associated to these corrections are at the sub- 
percent level relative to the total antineutrino spectrum. 
C. Impact on published measurements 

The new conversion of the ILL data described in sec- 
tion IIVI leads to a 2.5% increase of the detected an- 
tineutrino flux while the predicted shape is basically un- 
changed. The impact of this correction on the analysis 
of published measurements of antineutrino oscillations at 
reactor is discussed by the authors in a separate arti- 
cle [HHl- The sensitivity of forthcoming reactor experi- 
ments is also updated in this context. 



VI. CONCLUSION 



FIG. 13. Variations of tlie U antineutrino spectrum for dif- 
ferent irradiation times with respect to a reference spectrum 
considered at equilibrium. 



Note that our reference spectrum at equilibrium does 
not use cumulative yields. Instead it is computed by 
our evolution code using independent fission yields and 
a long irradiation time. In the analysis of cumulative 
yields it is in fact assumed that all nuclei have reached 
equilibrium. This is fully justified for short lived fission 
products, whereas there is some ambiguity with the decay 
products of longest half-lives. To avoid apparent double 
counting of the cumulative yields of the daughters, some 
very long-lived decays have been removed from the cumu- 
lative yield databases but some are still present in the li- 
braries. Because of these problems, it is recommended by 
nuclear databases (53| as a safer method to use indepen- 
dent yields with an inventory code. We have considered 
450 days of irradiation as our reference to account for a 
spectrum that would have reached quasi-equilibrium. In 
addition, this duration corresponds to a typical irradia- 
tion time of a fuel assembly in a PWR core. 

Off-equilibrium effects where also computed in (53 | 
where the authors used beta branches of 571 fission 
fragments. The fission yield were taken from j55| and 
beta decay properties came from experimental data. 
Our results are compatible with theirs considering the 
quoted uncertainties and the possible discrepancies in 
the neutron energy spectrum and flux used in the cal- 
culations. Note also that small additional discrepancies 
could arise from the smaller number of fission products 
used in the calculation of Kopeikin et al. 

In conclusion, our new reference spectra presented in 
section lTVl are. strictly speaking, valid only for irradiation 
times comparable to the ones used at ILL for their mea- 
surement. For longer irradiations, corrections to these 
spectra are listed in table I VIII The above comparison 
between independent estimates suggest that the system- 



Using all available data on fission yields and beta 
decays of fission products we have shown that ab initio 
calculations of total beta spectra agreed with the 
reference ILL beta spectra at the 10 % level, illustrating 
the tremendous amount of nuclear data collected today. 
From this work we gave a prediction of the antineutrino 
spectrum associated with the fission of ^ssy ^[f^^y 
estimated relative uncertainty increasing from 10 to 
20% with energy in the 2-8 MeV range. Since this 
isotope contributes to about 10% of the total fission 
rate of a reactor such a prediction is valuable. However, 
for the dominant isotopes, the remaining systematic 
errors of nuclear databases as well as the contribution 
of poorly known beta transitions still prevent the ab 
initio approach from any use for high precision neutrino 
oscillation experiments at reactors. 

This motivated the development of a new mixed- 
approach combining the accurate reference of the ILL 
electron spectra with the physical distribution of beta 
branches provided by the nuclear databases. We pre- 
sented how this method gets rid of the main systematic 
error of the ab initio approach allowing a better control 
of the conversion of reference electron spectra into 
antineutrino spectra. While the final error budget ended 
up being comparable to previous reference work [lol - [l^ , 
we demonstrated that the antineutrino spectra emitted 
by the fission of ^"^^U, ^^^Pu and ^"'^Pu isotopes have 
to be corrected for a systematic shift of about 3% in 
normalization. This net effect was presented as the 
combination of an improved implementation of finite 
size corrections to the Fermi theory plus a more realistic 
description of the distribution of /^-branches. 
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